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Abstract 

We study the quantum summation (QS) algorithm of Brassard, H0yer, Mosca and 
Tapp, see £Q, that approximates the arithmetic mean of a Boolean function defined 
on N elements. We improve error bounds presented in in the worst-probabilistic 
setting, and present new error bounds in the average-probabilistic setting. 

In particular, in the worst-probabilistic setting, we prove that the error of the QS 
algorithm using M — 1 quantum queries is ^ttM^ 1 with probability which improves 
the error bound 7rM _1 + tt 2 M~ 2 of pQ. We also present error bounds with probabilities 
p G (|, and show that they are sharp for large M and iVM -1 . 

In the average-probabilistic setting, we prove that the QS algorithm has error of 
order min{M _1 , iV -1 / 2 } iff M is divisible by 4. This bound is optimal, as recently 
shown in JHJ- For M not divisible by 4, the QS algorithm is far from being optimal if 
M <C N 1 / 2 since its error is proportional to M _1 . 

1 Introduction 

The quantum summation (QS) algorithm (also known as the amplitude estimation algorithm) 
of Brassard, H0yer, Mosca and Tapp, see pQ , computes an approximation to the arithmetic 

*This research was supported in part by the National Science Foundation (NSF) and by the Defence 
Advanced Research Agency (DARPA) and Air Force Research Laboratory under agreement F30602-01-2- 
0523. 
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mean of all values of a Boolean function defined on a set of iV = 2 n elements. Information 
regarding the Boolean function is supplied by quantum queries. The quantum queries play a 
role similar to the use of function values in the worst case and randomized settings. Suppose 
that we use M — 1 quantum queries. Obviously, the only case of interest is when M is much 
smaller than N. It was proven in PQ that the error of the QS algorithm is at most 

Jf + Jp with probability ^ = 0.81 ... . (1) 

Nayak and Wu, see [7j, showed that for any p £ (|, 1] the error of any quantum algorithm 
that uses no more than M — l quantum queries must be proportional to M _1 with probability 
p. Therefore, the QS algorithm enjoys the smallest possible error modulo a factor multiplying 

The minimal error estimate of order M~ l in the quantum setting should be compared to 
the minimal error estimates in the worst case and randomized settings of algorithms using 
M — l function values. It is known, see [S], that in the worst case setting, the error bound is 
roughly |(1 — M/N). This means that as long as M is much less than N the error is almost 
|, and is therefore of order M times larger than in the quantum setting. In the randomized 
setting, the classical Monte Carlo is almost optimal, and the error bound is roughly 1/(2 vM), 
see again jSj. Hence, it is of order \[M larger than in the quantum setting. 

The QS algorithm has many applications. In particular, it can be used for approximation 
of the arithmetic mean of a real function, which is the basic step for approximation of many 
continuous problems such as multivariate integration, multivariate approximation and path 
integration, see 01 HE EI] . 

Since the QS algorithm has so many applications, it seems reasonable to check whether 
the estimate (JTJ is sharp and how the error decreases if we lower the probability p = to 
p > |. It also seems reasonable to study the error of the QS algorithm in various settings. 
The estimate (0) corresponds to the worst-probabilistic setting, which is most frequently 
used in the quantum setting. The essence of this setting is that it holds for all Boolean 
functions. It is also interesting to study the average performance of the QS algorithm with 
respect to some measure on Boolean functions. This is the average-probabilistic setting. In 
the worst-average and average-average settings, we study the worst or average performance 
with respect to Boolean functions and the average performance with respect to all outcomes 
of a quantum algorithm. We add in passing that the worst-average setting is usually used 
for the study of the classical Monte Carlo algorithm. 

Sharp error bounds in the worst- and average-probabilistic settings are addressed in this 
paper whereas the worst- and average- average settings will be studied in a future paper. We 
study error bounds with probabilities p £ (~, -%]■ If we want to obtain error bounds with 
higher probability, it is known that it is enough to run the QS algorithm several times and 
take the median as the final result, see e.g., jlj. 

In the worst-probabilistic setting, we show that ((TJ can be slightly improved. Namely, 
the error of the QS algorithm is at most 
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Furthermore, for large M and N/M we prove that the last estimate is sharp. More 
generally, for p e (~, we prove that the error of the QS algorithm is at most 

(1 - v- l {p))-n . , _ i ... , 

with probability p, 



M 

where w -1 is the inverse of the function v (A) = sin 2 (7rA)/(7rA) 2 . We prove that the last esti- 
mate is sharp for large M and N/M. We have 1 — e (|, |] and it is well approximated 
by jq^ 2 P + i- In particular, for the most commonly used values of p we have 



v 



- 1 (|+))vr = 1.75 . . . , (1 - v-\l))n = 2.23 . . . , (1 - v -\$)) = f tt = 2.35 .... 



In the average-probabilistic setting, we consider two measures on the set of Boolean 
functions. The first measure is uniform on Boolean functions, while the second measure is 
uniform on arithmetic means of Boolean functions. The results for these two measures are 
quite different. The mean element of the arithmetic means is ~ for both measures. However, 
the first moment is of order iV -1 / 2 for the first measure, and about | for the second. The 
first moment is exactly equal to the error of the constant algorithm that always outputs |. 
This explains why we can obtain the error of order A r ~ 1 / 2 without any quantum queries for 
the first measure. This provides the motivation for us to check whether the error of the QS 
algorithm enjoys a similar property. It turns out that this is indeed the case iff M is divisible 
by 4. That is, for M divisible by 4, the average-probabilistic error of the QS algorithm is of 
order min{M _1 , iV" 1 / 2 }, and if M is not divisible by 4, then the error is of order M _1 . For 
the second measure, since the first moment is not small, the average-probabilistic error of the 
QS algorithm is of order M _1 for all M. For both measures, the upper bounds presented in 
this paper match lower bounds that were recently obtained by Papageorgiou, see fjjj . Hence, 
the QS algorithm enjoys minimal error bounds also in the average-probabilistic setting if we 
choose M divisible by 4 for the first measure and with no restriction on M for the second 
measure. 

The quantum setting, and in particular the QS algorithm, is relatively new and probably 
not well known, especially for people interested in continuous complexity. Hence we present 
all details of this algorithm, emphasizing its quantum parts. Since we wanted also to find 
sharp error bounds, we needed a very detailed analysis of the outcome probabilities of the 
QS algorithm. 

We outline the contents of this paper. In Section El we define the QS algorithm. Section 
El deals with the performance analysis of the QS algorithm in the worst-probabilistic setting, 
see Section 13 .l\ and in the average-probabilistic setting, see Section 13.21 



2 Quantum Summation Algorithm 

We consider the most basic form of the summation problem, i.e., the summation of Boolean 
functions. Let M N denote the set of Boolean functions / : {0, . . . , N — 1} — > {0, 1}. Let 

a / = iv E /(*) 

i=0 
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denote the arithmetic mean of all values of /. Clearly, aj G [0,1]. 



Problem: For / G M N , compute an ^-approximation aj of the sum aj such that 

\(if — cif \ < e. (2) 

We are interested in the minimal number of evaluations of the function / that are needed 
to compute df satisfying (j2j). It is known that in the worst case setting, we need roughly N(l— 
e) evaluations of the function /. In the randomized setting, we assume that df is a random 
variable and require that (J2J) holds for the expected value of \df — o/| for any function /. It is 
known, see e.g., [Hj, that in the randomized setting we need roughly mm{N, e^ 1 ^ 2 } function 
evaluations. In the quantum setting, we want to compute a random variable df such that 
(J2J) holds with a high probability (greater than |) either for all Boolean functions or on the 
average with respect to a probability measure defined on the set Bjy . These two error criteria 
in the quantum setting will be precisely defined in Section 01 

In this section we describe the quantum summation algorithm, which is also called the 
quantum amplitude estimation algorithm. This algorithm was discovered by Brassard, H0yer, 
Mosca and Tapp [T], and uses Grover's iterate operator as its basic component, see j2J. We 
use standard notation of quantum computation, see e.g., 

For simplicity we assume that N = 2™. Let Ti n denote the tensor product C 2 (g> ■ ■ ■ <g> C 2 
of n copies of C 2 , with C 2 the 2-dimensional complex vector space. Unit vectors from C 2 are 
called one qubit quantum states (or qubits). Let |0) and |1) be an orthonormal basis of C 2 . 
Then any qubit can be represented as 

|V) = '0o|O)+Vi|l) with and |V> | 2 + hM 2 = 1- 

For j = 0, 1, . . . , N - 1, we have j = YT k Zl ^^Jk, with j k G {0, 1}. Let 

n-l 

\j) = ®\h). 

The set {\j) : j = 0, . . .,N — 1} forms an orthonormal basis of Ti, n and any unit vector 
\ip) G 7i n can be represented as 

N-l N-l 

\^) = ^2^j\j) with and ^|^f = l. 

3=0 3=0 

Unit vectors from TC n are called n qubit quantum states (or quantum states or just states, 
whenever n is clear from the context). 

The only transforms that can be performed on quantum states are defined by certain 
unitary operators on Ti n . We now define the six unitary operators that are basic components 
of the summation algorithm. Since unitary operators are linear, it is enough to define them 
on the basis states 

1. Let Sq : Tin Ti n denote the inversion about zero transform 

So|j> = (-l)^|j>, 
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where 5j t0 is the Kronecker delta. Hence, S |0) = — 10) and Sob') = b) f° r a ^ 3 7^ 0- 
This corresponds to the diagonal matrix with one element equal to —1, and the rest 
equal to 1. The operator So can be also written as the Householder operator 

S = /-2|0)(0|. 

Here, for a state we let denote the projection onto the space span{|^)} 

given by 

(iv>)<v>i)k) = <#zm 

where (ip\x) is the inner product 1 in TC n , (ip\x) = V^fc- The matrix form of the 

projector |'0)('0| m the basis {\j)} is (fpk^^f^o- O ne can a l so view the matrix form 
of the projector \ip){ifj\ as the matrix product of the Af x 1 column vector and the 
JV x 1 row vector which is the Hermitian conjugate of = \^)- For any 

\x) G 7i n we have 



|(/-2|0)(0|)|x) = (k\x) — 2<0|aj><A:|0> 



Xk — 2xfc = — Xk for k — 0, 
Xk — = Xk for k 7^ 0. 



Hence / — 2|0)(0| = So, as claimed. 
2. Let Wn '■ 7~tn — > 7~L n denote the Walsh- Hadamard transform 

n—l 

^b') = ^=(g)(|0) + (-1)^|1)). 

^ N fc=0 

That is, the Walsh- Hadamard transform corresponds to the matrix with entries 

(i\w N \j) = ^=li(i k \ (|o) + (-i)^|i)) = -L Jj(-i)** = -L(-i)^ 1 *"*. 



\N-1 



The matrix ((i\WN\j))i J=0 is symmetric. Furthermore, 



•,j= 

n-l 



w N \j) = J=w B <g)(|o> + (-i)*|i)) 

1 'X 1 / 1 r_i v'fc 

= ^®b (|0)+,1))+L i- ,|0> - |1)) 

Thus, W 7 "^ = / and W^ 1 = Wn is orthogonal. This means that the operator Wn is 
symmetric and unitary. 



1 We follow the quantum mechanics notation in which the first argument is conjugated in the inner product, 
whereas in the standard mathematical notation the second argument is usually conjugated. 
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F K ,n\j) 



3. For K = 1, 2, . . . , 2 n , let F Kn : H n — > 7i„ denote the quantum Fourier transform 

T^k=o^ ijk/K \k), for j=0,l,...,^-l, = 
|j> for ./ A' 2"- 1. 

Hence, F Kn corresponds to the unitary block-diagonal matrix 

F K 
/ 

where Fk = (K ~ x I 2 e 2m ^ k ^ k ^ k ^ is the matrix of the inverse quantum Fourier trans- 
form. For K = 2 n = N we have 



N-l 1 N-l /N-l \ 

7 nM = Y,^F N;n \j) = iY,^ e2m3k/N ) i*>- 



j=0 v k=0 \j=0 

The coefficients of Fjv, n |'0) m the basis {|j)} are the inverse quantum Fourier transforms 
of the coefficients of the state Note that Wn and F NjH coincide for the state |0), 
i.e., 

N-l 

w N \o) = F N>n \o) = -7=5]b'). 



N ■ n 
3=0 

4. Let Sf : 7i n — > 7Y n denote the quantum query operator 

s f \j) = (-i) m \j), 

This again corresponds to the diagonal matrix with elements ±1 depending on the 
values of the Boolean function /. This operator is the only one that provides informa- 
tion about the Boolean function /. This is analogous to the concept of an oracle or a 
black-box which is used in classical computation and which supplies information about 
the function / through its values. 

The standard definition of the quantum query Sf is 

s f : n n ® c 2 - n n ® c 2 , s f \j) \i) = \j) \i © f(j)), 

where © means addition modulo 2. We can simulate Sf by Sf if we use an auxiliary 
qubit (l/v / 2)(|l> - |0)), namely, 



U) li)-io) ^_ U) lie/(i))-i/Q-)) 



V2~ J V? 



( _ 1)m|)) o>. (S;lj)) oi. 



5. Let Qf : H n — > H n denote the Grover operator 

Q f = -W N S W N 1 S f . 
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This is the basic component of Grover's search algorithm, see |2|. As we shall see, 
Qf also plays the major role for the summation algorithm. The eigenvectors and 
eigenvalues of Qf will be useful in further considerations. Let 



1 N ^ 



and j^o), denote the orthogonal projections of \ip) onto the subspaces 

span{|j) : f(j) = 0} and span{|j) : f(j) = 1}, respectively. That is, 



AT 

k: f(k)=j 

Then |-0) = \ip ) + \ipx) and (^ol^i) = 0. Furthermore, (ipj\ipj) = N' 1 J2k-.f(k)=j for 
j = 0, 1, so that (ipi\ipi) = a and (00 1 "00 ) = 1 — a, where a = a/ is the sum we want to 
approximate. 

From pQ, we know that 

Q / |^o) = (l-2a)|^o) + 2(l-a)|^ 1 ), 

Q/IV-i) = -2a|0 o ) + (l-2a)|^i). 1 J 

For the sake of completeness, we provide a short proof of J3J). By the definition of the 
operator we have 

s f \^) = (-iy\^), j = o,i, 

and 

W^SoW^ 1 = W N (I-2\0)(0\)W„ 1 = 7-2(W^|0)(0|W^). 
Since (0\W N = (W N \0)Y = (|^)) t = (ip\, we obtain for j = 0, 1, 

WsSoW^j) = |^} -2(|V)(Vi) |^> = hfc) -2(^)1^) = |Vi> -2(Vil^)l^>- 
From this we calculate for j = 0,1, 

<2j#i) = (-i) 1+ ^s w^V;> 

= (-l) 5 ^ (|^> - 2(* ifl a + S jfi (l - a))(|^o) + |V>i») , 

which is equivalent to (j3J). 

Thus, the space span{|-0o), |^i)} is an invariant space of Qf and its eigenvectors and 
corresponding eigenvalues can be computed by solving the eigenproblem for the 2x2 
matrix 

1 - 2a -2a 
2(1 -a) I -2a 

For a G (0, 1), the eigenvectors and the corresponding orthonormalized eigenvalues of 
Qf are 

|^ ± ) = ^r±-^L=|^ > + ^|V'i>J and \ ± = l-2a±2i^a(l-a)=e ±2i9 % 

7 



where 9 a = arcsin y/a. Moreover, it is easy to check that 

m = ^(e i6 ^ + )-e- ie ^)). (4) 

For a G {0, 1}, we have span{|^ ), \ipi)} = span{|^)} and \ip) is the eigenvector of Qf 
with eigenvalues ±1, respectively. For a G {0, 1}, we define 

= i X ~ a V2\il)) and = 0. 

Then it is easy to check that (jU is valid, and \± = e ±2tda = (— l) a is an eigenvalue of 
Q ; for all a G [0,1]. 

6. The next unitary transform, called the Grover iterate operator, is defined on the tensor 
product of H m <%> 'Hn and uses m + n qubits. The first space Ti m and m qubits will 
be related to the accuracy of the quantum summation algorithm. The Grover iterate 
operator A m (Qf) : 7i m ® 7i n — > Ti. m ® Ti n is defined by 

Kn(Q f )\j)\y) = \j)Qj\y) for \j)\y) e H m ® H n . 

Hence, the power of Qf applied to the second component depends on the first one. 
Note that j may vary from to 2 m — 1. Therefore A m (Qf) may use the powers of Qf 
up to the (2 m - l)st. 



We need one more concept of quantum computation, that of measurement. Suppose s is 
a positive integer and consider the space TL S . Given the state 



X^fclfc) G H s 



k=0 

we cannot, in general, recover all the coefficients ipk- We can only measure the state with 
respect to a finite collection of linear operators {Mj}? =0 , where the Mj : 7i s — > 7i s satisfy 
the completeness relation 

3=0 

After performing the measurement, we obtain the outcome j and the state collapses into 
the state 



1 



M 3 \ 



'{^\M)M^) 

both these events occur with probability (ip\Mj Mj\ip) . Note that for Mj\if>) = the outcome j 
cannot happen with positive probability. Hence, with probability 1 the outcome j corresponds 
to Afj-|^) ^0 for j = 0,l,...,p. 

The most important example of such a collection of operators is {b')0'|},-="o ■ Then, 
the measurement of the state with respect to this collection of operators gives us the 
outcome j and the collapse into the state 
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with probability \ipj\ 2 , j = 0, 1, . . . , 2 s — 1. 

Another example is a variation of the previous example and will be used in the quantum 
summation algorithm. We now let s = m + n, as for the Grover iterate operator, and define 
Mj :H m ®H n ^ H m ® H n by 

for j — 0, 1, . . . , 2 m — 1, with I denoting the identity operator on H n . That is, 

(b'}0'|® J) \x)\y) = (j\x)\j)\y) 

for \x) G H m and \y) G H n . 

Since X^o" 1 (U)O'I ® -0 \ x )\v) = \ x )\v) f° r a ^ basis states \x) of Tt m and |y) of 7i n , the 
completeness relation is satisfied. Consider now the probability of the outcome j for a special 
state of the form \ip) = l^i)!^) with l^i) e H m , \ip2) G 7^„ and (^fcl^fc) = 1 for = 1,2. 
Since <S> I is self-adjoint, the outcome j and the collapse of the state \ip) to the state 



occur with probability |(j|^i)| 2 . Hence, this collection of operators measures the components 
of the so-called first register 1^) of the quantum state \ip). 

Following |Tj. we are ready to describe the quantum summation (QS) algorithm for solving 
our problem. The QS algorithm depends on a Boolean function / and on an integer parameter 
M that controls the number of quantum queries of / used by the algorithm. We perform 
computations in the space 7i m <8> TC n , with m = |~log 2 M] , so we use n + m qubits. As we will 
see later, the accuracy of the algorithm is related to the dimension of the space Ti. m . 

Algorithm QS(/, M) 

Input state: |0)|0) G H m ® H n with m = [log 2 M] and n = log 2 N. 
Computation: 

1. \Vi) = F M , m ®W N \0)\0), 

2. \v2)=A m (Q f )\rn), 

3- \v 3 ) = (F^ m ®I)\ m ). 



Measurement : 

Perform the measurement of the state [773) with respect to the collection 



{ (li)OI) ® I }?=o 1 - Denote the outcome by j. 



Output: af(j) = sin 2 (irj/M). 



We briefly comment on the QS algorithm. The input state is always the same and does 
not depend on /. Step 1 computes ]%) = {NM)^ 1 / 2 Yl!j^ Y^kZo which is the equally 

weighted superposition of the basis states. Step 2 computes |^ 2 ) by using the Grover iterate 



9 



operator. During this step we use the successive powers of the Grover operator Qf, and this is 
the only step where information about the Boolean function / is used. We shall see that the 
QS algorithm uses M — 1 quantum queries. Step 3 computes by performing the inverse 
quantum Fourier transform on the first m qubits, and prepares the system for measurement. 
After Step 3, we perform the measurement, obtain the outcome j and compute the output 
a,f(j) on a classical computer. We stress that the distribution of the outcomes j depends on 
the Boolean function /, and this is the only dependence of the output a,f(j) on /. 

3 Performance Analysis 

In this section we analyze the error of the QS algorithm. As we have seen in Sectional the 
output df(j) of the QS algorithm is a random value chosen according to a certain distri- 
bution dependent on the input function /. In this way, the QS algorithm is a randomized 
algorithm. Various ways of measuring the performance of randomized algorithms are com- 
monly used in the analysis of algorithms and computational complexity. They correspond 
to various error criteria. In this paper we consider two error criteria: worst-probabilistic and 
average-probabilistic. In a future paper we consider other two error criteria: worst-average 
and average- average, which correspond to the worst or average performance with respect to 
Boolean functions and the average performance with respect to all outcomes. 

Worst-Probabilistic Error 

We start with the error criterion that is used in most papers dealing with quantum computa- 
tions. We are interested in the worst case error of the QS algorithm that holds with a given 
probability p. Here p G [0, 1] and 1 — p measures the probability of QS algorithm's failure and 
usually p is set to be |. In our analysis, however, we will allow an arbitrary p G (|, -%]. The 
choice of the upper bound = 0.81 . . . will be clear from the analysis of the QS algorithm. 
The QS algorithm outputs df(j) with probability Pf(j) for j = 0, 1, . . . , M — 1, see Theorem 
121 where the Pj(/)'s are given. Its worst-probabilistic error is formally defined as the smallest 
error bound that holds for all Boolean function with probability at least p, i.e., 



j: \a f -a f (j)\<a 

It is easy to see that e wor_pro (M, p) can be rewritten as follows. Let A C {0, 1, . . . , M — 1}. 
For / G Btv define the measure of A as 



e 



, wor— pro 




MA/) 



Then 



e 



wor— pro 



(M,p) 



= max 

f&N A: 




(5) 
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Average-Probabilistic Error 



The worst-probabilistic error e wor ~ pro (M, p) of the QS algorithm is defined by the worst 
performance with respect to Boolean functions. It is also natural to consider the average 
performance of the QS algorithm with respect to Boolean functions. Let p be a probability 
measure on the set B^v- That is, any Boolean function / G B^v occurs with probability p(/). 
Obviously, p(/) > and J^ eBjv p(/) — 1- The average-probabilistic error is defined by 
replacing the first max in (jSJ) by the expectation, i.e., 

e avg- P ro (M)p) = J2 p(/) min msx\a f - a f (j)\, 

^— ' A: n{AJ)>p J&A 



Hence, we are interested in the average error that holds with a certain fixed probability. 



3.1 Worst-Probabilistic Error 

We begin by citing a theorem from pQ for which we will propose a number of improvements. 

Theorem 1. rij/ For any Boolean function f G B^, the QS algorithm uses exactly M — 1 
quantum queries and outputs a that approximates a — a/ such that 

, . 2-K /—. 7T 2 7T 7T 2 

a — a < — vafl — a) H r < 1 ~ 

with probability at least A = 0.81 .... Hence, 

E p/0') > 4 v ^ G 

Z — c ^ 

j: |a / (j)-a / |<(27r/M) v /a / (l-a / )+7r2/A/ 2 

and, therefore, 



7T 7T 2 



v '- 2; -m M 2 

Using proof ideas of Theorem ^ from pQ we present the following theorem and the subse- 
quent corollaries. 

Theorem 2. For any Boolean function f G M N , denote 

o~ a = Pat — — arcsin Ja G [0, \M\ . 

1 7T 

1. TTie algorithm uses exactly M — 1 quantum queries, and log 2 iV + |~log 2 M] qubits. 

2. For j — 0, 1, . . . , M — 1, the outcome j of the QS algorithm occurs with probability 

S m 2 (n(j-a af )) ( sin 2 (tt(j - cr af )/M) \ 
Pf[3) 2M 2 sin 2 ^' -a af )){ + sin 2 (n(j + a a ,)/M) J ' W 

(If sin(7r(j ± <7 a ,)/M) = we need to apply the limiting value of the formula above.) 
For j = M, M + 1, . . . , 2^ log2 M l — 1 ; £/ie outcome j occurs with probability 0. 
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3. If a af is an integer then the QS algorithm outputs the exact value of aj with probabil- 
ity 1. This holds iff a/ = sin 2 (/c7r/M) for some integer k G [0, \M\. In particular, this 
holds for af = 0, for a/ = 1 and even M , and for a/ = \ and M divisible by 4. 

4- Let x = 7r(|~<7 a ] — <7 a )/M and x = n(a a — [a a \)/M. If a af is not an integer then the 
QS algorithm outputs the same value a = a/(|~a a ]) = a/(M — |~cr a ~|) for the outcomes 
\a a ~\ and M — \a a ~\ such that 



a — a = 



sin(:r) ^2 a/ a(l — a) cos(:r) + (1 — 2a) sin(:r)j 



(7) 



with probability 



sin 2 (7r([o- a ] - a a )) 
MW(£(K| - a a )) 



sin 



, . j f , .J(7r(ra Q l-a a )/M) ^ 

sm (vr([a a l + a a )/M)J 



sin 



> ^(^,^2 =l- Y (M-*a) +0((\a a ]-a a )), 



(8) 



and outputs the same value a = a/(|_°"aj) = — ^L CT aJ,o)-^ — / or ^ e outcomes 
|_o" a J <J^<^ (1 — ^L°"aJ:o)^ — [cr a \such that 



a — a 



sin(x) ^2a/o(1 — a) cos(rr) + (1 — 2a) sin(V)j 



7T 



< ^K-L^aJ) 



(9) 



with probability 



sin 2 (7r(o- a - \a a \)) 
M2sin 2 (^K-La a J)) 



l + (l-^KJ,o) 



sin 2 (7r( 



a a -[a a \)/M) \ 



sm 2 {7i{a a +[a a \)/M)J 



> 



Sill 7T <7 



7T 2 ((T a - [OaJ) 



^# = i - £ K - M) 2 + o(K - La a j) 4 ). 



(10) 



Proof. As before, let # a = arcsin y/a and 

M-l 
fc=0 

for arbitrary w6K. Note that 



1 M-l 

|5 M (c)) = ^=^ e 2 - fe |A;), i = </=!, 



\S M (j/M)) for; = 0,1,..., M-l, 

|j) for j = M,M + l,...,2 m -l. 
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The steps 1-3 of the QS algorithm are equivalent to the application of the operator (F Mm <g> 

I) A m (Q f) FM,m ® Wn to the state |0)|0) G 7Y m <8> Ti n . Then |r/i) can be written as 
M _i /2 fM-i |'^|^^ and |^ = ig giyen by Rence 



l^i) 



—i 



M-l 



'2M 



£|j> (e^^ + )-e-^|^-»- 



j=0 



Applying A m (Q/) in Step 2 and remembering that QjlV'i) = we obtain 

M-l 



|»/2> = Am(Q/) 



—2 



'2M 



"2ijda e -id a 



l^->) 



j=0 



—2 

71 



(e*-|SW(0- a /M))|V + ) - e-^I^C-a-a/M))^-)) 



Since j = 0, 1, . . . , M — 1, the largest power of Q/ is M — 1. Hence, we use exactly 
M — 1 quantum queries to compute |r] 2 ). The remaining steps of the QS algorithm do not 
use quantum queries. This means that the total number of quantum queries used by the QS 
algorithm is M — 1, and obviously we are using n + m qubits. This proves the first part of 
Theorem El 

Step 3 yields the state 



\Vs) = (F M ] m ®I)\ m ) = -± {^F M ]JS M {a a /M))\i, + ) - e^F^JS M {-a a /M))\^-)) ■ 

We are ready to analyze the probability of the outcome j of the QS algorithm. Observe that 

\a±) ■= F M ] m \S M {±a a /M))\^ ± ) 
= (j\F M ] m \S M (±a a /M))\j}\^ ± } 

(S M (j/M)\S M (±a a /M)) \j)\ij ± ) for j = 0, 1, . . . , M - 1, 
for j = M,M + l,...,2 m -l, 



and therefore 

(a±|a±) 



\(S M (j/M)\S M (±a a /M))\ 2 (^±|V±> for j = 0, 1, . . . , M — 1. 







for j = M,M + l,...,2 m -l. 



Observe that for a G (0,1), we have (i[>±\if)±) = 1, whereas for a G {0,1}, we have 
<V+IV>+)= 2and(^_|^) = 0. 
For w 1 ,w 2 6 M we have 



|<SmM|SmM>| : 



1 

M2 



1 M-l 

M 

3=0 
M-l 

3=0 



M-l 



01 



M ^ 

3=0 



-27ri(£Ji— 0*2)3 
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If u\ — lu 2 is an integer then the last sum is clearly M, and the whole expression is 1. If 
uj\ — uo 2 is not an integer then 



I ^Z} , . -27riM( Wl -o;2) _ 1 

_ p -27ri(o; 1 -^ 2 )j r 1 

m 2^ e 



M j-^ M (e- 2 ^(^- W2 ) - 1)' 

which holds for all loi,u 2 G M if we take 0/0 as 1. Therefore 



M-l 



j=0 



1 — cos(27rM(c<Ji — lu 2 )) sin (vrM^ — uo 2 )) 



M 2 (l - cos(2tt(u;i - u 2 ))) M 2 (sin 2 (7r(^i - w 2 ))) 



Hence 



2 sin 2 (M7r(co'i — CJ2)) 
M 2 sin 2 (7r(wi - u 2 )) 



which holds for all u 2 G K. if we take 0/0 as 1. Applying this we conclude that 

(a± i Q±) = f «.a%&:W **fe> for < " ' * 

[0 for j = M,M + l,...,2 m - 1. 

The outcome j occurs after the measurement and the state |r/ 3 ) collapses to the state 
((%|mJm,-|77 3 ))- 1 M,-|77 3 ), where Mj = <g> I. For j = 0, 1, . . . , M - 1, we have 



e -i0a 



(S M (j/M)\S M (-a a /M))\^)) 



whereas Mj|^ 3 ) = for j = M, M + 1, . . . , 2 m — 1. Since and \ip-) are orthogonal we 
have 

rt 



(TftlMjM^s) = |((a+K> + («-!«-))■ 
Hence, the outcome j, j = 0, 1, . . . , M — 1, occurs with probability 

/ ■ \ = I ( sin 2 (7r(j - a a )) sin 2 (7r(j + a a )) \ 

PfU) 2\M 2 sm 2 (7i(j -a a )/M) M 2 sm 2 (7i(j + a a )/M)J ' 1 ' 

Indeed, for a G (0, 1), we have (V'ilV'i) = 1 an d follows from the form of {oi±\a±). For 
a G {0, 1}, we have = 2 and (^-l^-) = 0. Since the two terms in (|12jl are now the 

same, the formula for (a + \a + ) again yields (fT2"J). 

Since sin 2 (7r(j — cr a )) = sin 2 (7r(j + cr a )), the last formula is equivalent to (0). Obviously 
for j = M, M + l, . . . ,2 m — 1, the probability of the outcome j is zero. This proves the second 
part of Theorem |21 

Assume now that a a G Z. If <r a = or er a = |M (if M is even) then the probability 
Pf{c a ) of the outcome a a is 1. For o a = we have a = and the output is a/(0) = 0. 
For cr a = 7;M we have a = 1 and the output is o/(|M) = 1. Hence, in both cases the QS 
algorithm outputs the exact value with probability 1. 
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If a a G Z and a a ^ {0, |M} then the probability of the distinct outcomes a a and M — a a 
is ~. These two values of the outcomes yield the same output 

sin 2 (vra a /M) = sin 2 (tt(M - a a )/M) = a. 

Hence, the QS algorithm outputs the exact value with probability 1. This proves the third 
part of Theorem El 

We now turn to the case when a a ^ Z. It is easy to check that the third part of Theorem [21 
holds for M — 1. Assume then that M > 2 which implies that \\M\ < M — 1. Since a a is 
not an integer, we have |~er a ] > 1, |~er a ] < | ^Mj < M — 1 and M— |~<x a ] < M — 1. This means 
that both |~cr a ] and M — \a a ] may be the outcomes of the QS algorithm. Obviously, these 
two outcomes are different iff \a a ] ^ |Af. Similarly, both [a a \ and (1 — S^j^M — [a a \ 
may be also the outcomes. They are different iff \_cr a \ ^ 0. 

We show that the outputs for the outcomes \a a ] and [a a \ satisfy (0) and Q with prob- 
abilities (jHJ) and ([EH) , respectively. We focus on the output for the outcome [cx a ] and its 
probability. The proof for the outcome [a a \ is similar. 

We estimate the error of the QS algorithm for the output a = sin 2 (7r [cr a ] /M). Recall 
that x = 7r( [cr a ] — a a )/M. We have 

|a - a| = | sin 2 (7r \a a ] /M) - sin 2 (7T(T a /M)| = | sin(x) sin(x + 2-Ka a /M)\ 
= | sin(x) ( sin(27rcr a /M) cos^) + cos(27ra a /M) sin(x)) | 
< 7r(rcr a l -a a )/M. 

Since sin(27ro" a /M) = 2a/o(1 — a) and cos(27rcx a /M) = 1 — 2a, this proves the estimate of the 
error of the QS algorithm in the fourth part of Theorem [21 

We find the probability of the output a. Since sin 2 (7rt/M) is injective for t G [0, |M], the 
output a occurs only for the outcomes |~cr a ] and M — \a a ~\ . If \a a ] = |M then these two 
outcomes are the same and a occurs with probability pj(^M). Due to (fT^) 

sin 2 (7r(iM-(T a/ )) 
PM2 ' M 2 sin 2 (vr(iM- ( T a/ )/M) 

which agrees with the claim in Theorem El 

If \a a ] ^ \M then \a a ] ^ M - \a a ] and a occurs for exactly two distinct outcomes. The 
probability of a is now equal to the sum of the probabilities P/([cr a l) +Pf(M — [cr a ] ) with 
P/'s given by (JT2J). Since both terms are equal, the probability of a is 2pf(\a a ]) which also 
agrees with the claim in Theorem Since sin(-^([cr a ] — a a )) < ^([cr a ] — a a ) we have 

sin 2 (7r([a a ] - a a )) sin 2 (7r \a a ] - a a ) 
M2sin 2 (f (\a a ] -a a )) ~ Tr*(\a a ] - a a )* ' 

We finish proving (jHJ) using the standard expansion of the sine. This completes the proof. □ 

Based on Theorem[2]we present simplified estimates of the error of the QS algorithm and 
of the corresponding probability. 
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Corollary 1. The QS algorithm outputs a such that 



IX 

\a ~ a\ < — m&x{\cr a ] - a a , a a - \a a \} (13) 



with probability at least 



sin (7r(r<7 ] - a a )) sin (7r(a a - [a a \)) 



.2/ vr /r_ 1 . \\ T )i,r2„:„2/ tt / _ i _ inn — _9 • 



AfW(£([Vr.l " *a)) M 2 sin 2 (f K - [a a \)) ~ ^ 

Proof. It is enough to prove Corollary ^ if a a is not an integer. The estimate of the error of 
the QS algorithm by the maximum of the estimates (J7J) and Q holds with probability that 
is the sum of the probabilities (jEJ) and |T0|) . Moreover, \a a ] — cr a = 1 — (a a — [o~ a \). It now 
suffices to observe that the function 

sin 2 (7rA) sin 2 (7r(l - A)) 
^ (A) = ^A^ + vr 2 (l-A) 2 

is a lower bound of the left hand side of ()14j) with A = \a a ] — a a , and attains the minimum 
^2 on the interval [0, 1] for A = |, see also PQ. □ 



Corollary Q guarantees high probability of the estimate (jlHjl . Unfortunately this estimate 
does not preserve the continuity of the estimates (J7J) and Q with respect to |~er a ] — cr a and 
fa — L°"aJ • The continuity of the estimates will be present in the next corollary at the expense 
of the probability of the outcome. This corollary will also play an essential role in the study 
of the average-probabilistic error of the QS algorithm. 



Corollary 2. The QS algorithm outputs a such that 



IX 

\a - a\ < —mm{\o- a ] - a a , a a - \a a \} (15) 



with probability at least 



f sin 2 (7r([a a l - a a )) sin 2 (7r(a a - [a a J)) 14 
maX \M 2 sin 2 (f ([a.] -a a ))'M 2 sin 2 (^K- KJ))J " vr 2 ' 1 > 

Proof. We may again assume that a a is not an integer. Let us define 

sin 2 (7rA) „ A 

Then u>(|~a a ] — a a ) is the probability of (J7j) and w(l — (|~er a ] — &a)) is the probability 
of P. For A G [0, §], note that w(-) is decreasing, and w(l — •) is increasing. Therefore 

w(A) > w{\) > w(l — A) for Ae[0,|]. 

Let |~er a ] — a a < a a — [a a \. Then \a a ] — a a < \. In this case (fi3j) is equivalent to (J7j) 
and holds with probability at least u>(|~o- a ] — cr a ), which corresponds to (fTBj) . Analogously, if 
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\cr a ] — a a > a a — [cr a \ then \cr a ] — a a > |. In this case (fT5|) is equivalent to © and holds 
with probability at least w(o~ a — [o~ a \), which also corresponds to Finally, note that 

max {w([<r a ] -a a ),w(a a - [a a \)} 

is minimal for \a a ~\ — cr a = ~ and is equal to 

1 . _ 2 7T 4 

sin ^* — 

M 2 2M - 7T 2 ' 

□ 



Unfortunately for |~cr a ] — <x a close to \ the probability of the estimate (fTo'j) is too small. 
However, in this case we may use Corollary^ which yields the estimate with high probability. 

We now turn to global error estimates, that is, estimates independent of a. Theorem^ of 
PP states, in particular, that \a — a\ < tt/M + tt 2 /M 2 with probability at least We now 
improve this estimate by combining the estimates (|13|) and (|15|) . 

Corollary 3. The QS algorithm outputs a such that 

la - a| < - — (17) 
AM K ' 



with probability at least That is, 



wor-pro/ M JM < £ 7T 

v 71 ~~ 4 M 



Proof. Let us define 



, x . sin 2 (7rA) sin 2 (7r(l - A)) ) 



Clearly, /i(|~cr a ] — cr a ) is a lower bound of the max{u>( |~<7 a ] — o~ a ),w(l — ([a ] — cr a ))} and 
therefore /i([cr a ] — cr a ) is a lower bound of the probability of the output satisfying (fTo'j) . We 
consider two cases. 

Assume first that A = \a a ~\ — a a G [0, |] U [|,1]. It is easy to see that then 
^i(A) > -^2 and the estimate (|15[1 yields 

\a - a\ < min{ \a a ] - a a , a a - [a a \ } <\jj 

with probability at least 

Assume now that |~er a ] — a a G (|, |). Then we can use the estimate JEJ), which holds 
unconditionally with probability at least 4j. In this case, we have 

|a - a| < max{ [a a ] - a a , a a - \a a \ } < - 
These two estimates combined together yield (fT7J). □ 
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The obvious consequence of Corollary |3] is that for M large enough we can compute the 
value of a exactly by rounding the output. 

Corollary 4. Assume that 

M>^N. 
2 

Then the rounding of the QS algorithm output to the nearest number of the form k/N yields 
the exact value of the sum a with probability at least 

The proof of Corollary |3] may suggest that the constant | in (|17j) can be decreased. 
Furthermore one may want to decrease the constant | at the expense of decreasing the 
probability These points are addressed in the next corollary. We shall see that the 
constant | may be lowered only by decreasing the probability. 



Corollary 5. Define 



and 



Then 



<7(p) = inf{c7: \a f -a f \<C^ V/ G 

sin 2 (7rA) 



D/V 



with probability at least p 



HA) = for A e [\, -]. 



for pe [0,4,), 
C{p)<{l-v-\p) for pe 

M/tt for pe (8/tt 2 ,!]. 



(19) 



Moreover, 1 — v 1 (p) G [|, |] and 

(l-^b)) 



2' 4J 
7T 2 1 

— V H — 
16^ 4 



< 0.0085 /or p G [ 



J_ JLl 

n-2 J ^2 



(20) 



Proof. For p G [0, -4-), Corollary El is a consequence of Corollary El Forp G 1], Corollary^] 
trivially holds since |o — a| < 1 = (M/tt)tt/M. For the remaining p's we use a proof technique 
similar to that of Corollary El 

Let p G [t 1 , -%]• It is eas y to check that v is decreasing and, therefore, v~ 1 (p) is well 
defined and w _1 (p) G [j, |]. We have to show that the estimate 



|a — a] < (1 - ^(p)) 



7T 

M 



(21) 



holds with probability at least p. We consider two cases. 

Assume first that A = \a a ] —a a G [0, v ~ 1 (p)} U [1 — v _1 (p), 1]. Observe that the function h 
defined in (118)1 can be rewritten as 



h(A) = max{w(A), u(l - A)}. 
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It is easy to see that in this case, h(A) > p and the estimate (fTo]) yields 

I 5 - a l < T7 min{A, 1 - A} < ^(p) -J- < (1 - ^(p))-^- 
i — M 1 J - M yF " M 

with probability at least p. 

Assume now that A = \a a ~\ — cr a G (v _1 (p), 1— ?; _1 (p)). Then we can use the estimate (|T3j). 
which holds unconditionally with probability at least ^ > p. In this case, we have 

|a - a| < ^ max{A, 1 - A} < (1 - ^ _1 (p))^- 

This proves (|2ip. 

We found the estimate (j^Oj) by numerical computations. □ 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 1: The estimate (fT^j) of C{p) for p G [0 



From Figure ^ we see that the estimate ()19j) is almost linear on the interval [4j, which 
explains why the right hand side of the estimate ()2U|) is small. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 2: The function v on [0, 1]. The two horizontal lines show -% and \ levels. The part 
of the graph between the arrows shows that v is almost linear. 



We now find a sharp bound on the worst-probabilistic error of the QS algorithm. We 
show that for large M and N/M the bound obtained in Corollary |S] is optimal for p e (|, -^]. 

Theorem 3. For large M and N/M, the worst-probabilistic error of the QS algorithm is 
given by 

e wor- P ro (M)p) = {1 _ v -i (p)) I_ {1 + 0{M -i ) + oiMN- 1 ) ) / orp e (±,£]. 
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Here, v is as in Corollary^ and 1 — v 1 (p) ~ (n 2 /16)p + | by \2(J\) . 

Proof. From Corollary it is enough to show a lower bound on the error. Define 



M j \ M 

For large M, we have 

Sl = \ + 0{M^) and 82 - 8l = sm(-)sm^ J 

There exist two Boolean functions /i and / 2 with sums ai = a/ x and a 2 = a/ 2 such that 

k - s< | < N~ x for z = l,2. 

Since c s . = \-gM\ + (i — 1) and the derivative of cx a for a = ^ is M/n, we have 

a ai = \\M] + (^(MiV^ 1 ) and a a2 = \\M] + 1 + O(MN^). 

Obviously, = ki/N for some integers ki with k\ < k 2 . Consider o~ x /n for x G {/ci, k\ + 
1, . . . , k 2 }. Then <t x /at varies from a ai for x = ki to <r a2 for x = k 2 . Since G [|, ~), 

for a positive and small rj (we finally let r\ go to zero), we can choose x = x v such that for 
a* = x v /N we have 

<r«. := [iMj + w - 1 (p)+r ? + 0(MiV- 1 ). 
For large N/M, we then have 

|/7 .J = |"±M] and = [±M] + 1, 

^a* - L^a*J = + 7? + 0(M/iV) and -a a .=l- v'^p) - V + 0(M/iV). 

Let a* denote the output for the outcome \a a *], and a 2 for Lc a *J- 
Due to (jZJ) and Q of Theorem we have 



a — ai 



a — a 



— * 



vf (1 - v-\p) - ri) (1 + 0(M _1 + MiV- 1 )) 
M (« _1 (p) + ry) (1 + C(M~ 1 + MN- 1 )) . 



Let us write 1 + o(l) for 1 + 0(r/ 2 + M _1 + MA^ 1 ). The probability of a* 2 is given by (PUJ 
and is now equal to 

sin^Tr^-ip) +rj)) , = sin 2 ^- 1 ^)) + ttt? sin(27rt>- 1 (p)) 

(7r(v _1 (p) + r?)) 2 1 1 JJ vr 2 t;- 1 (p) 2 (l + 2r//t;- 1 (p)) 1 1 JJ ' 

Since p = = sin 2 (7ru -1 (p))/(7n;""'' 1 (p)) 2 , the probability of a 2 is 

1 



/' ( 1 - ^1 ( " * cot^" 1 ^) } ) ( 1 i „( I ) ) 
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Since cot(t) < 1/t for t G [W, |7r], we see that the probability of is slightly less than p for 
small 77. 

We are ready to find a lower bound on the worst-probabilistic error 

e wor_pro (M, p) = max min max \at — a f f 7 ) I 

/eliv A:^(AJ)>p j&A 17 7WI 

of the QS algorithm. Take the function / that corresponds to a*. We claimed that the 
error is minimized if A = {[cr a *\ , [oa*]}- Indeed, [_a a *\ must belong to A since otherwise 
MA /) — 1 — P/(L°a*J) = 1 — P + o(l) < p for p > |. The probability of [cr a *\ is slightly 
less than p, and so the set A must also contain some other outcome j. If j = \a a *~\ then 
the error bound is roughly (1 — v^ 1 {p) — rj)ir/M, and the sum of the probabilities of the 
outputs for the outcome [a a *\ and \u a *] is always at least ■% > P- On the other hand, if 
\a a *1 does not belong to the set A then any other outcome j yields the output sin 2 (7rj/M). 
Since sin 2 (vr(j + 1)/M) - sin 2 (7rj/M) = sin(7r/M) sin(vr(2j + 1)/M), the distribution of the 
outcomes around | is a mesh with step size roughly tt/M. Hence, if j 7^ [o"a*l, the error is 
at least roughly (1 + v~ l {p))n /M > tt(1 — v~ 1 {p))/M. Thus the choice j = \a a *~\ minimizes 
the error and for rj tending to zero, the error is roughly (1 — v~ 1 (p))7r/M. This completes 
the proof. □ 

From these results, it is obvious how to guarantee that the error of the QS algorithm is 
at most e with probability at least p. Since |a — a] < (1 — v^ 1 {p))i[ /M holds with probability 
p, it is enough to take M > (1 — v ~ 1 (p))7r/e. Due to Theorem El this bound is sharp for small 
e and large eN. We have 

Corollary 6. For p e (|, the algorithm QS(f, [(1 — v^ 1 (p))n/e]) computes a with the 
error e and probability at least p with [(1 — v ~ 1 (p))n/e] — 1 quantum queries. For small e 
and large eN, the estimate of the number of quantum queries is sharp. 

3.2 Average-Probabilistic Error 

In this section we study the average performance of the QS algorithm with respect to some 
measure on the set B^v of all Boolean functions defined on the set {0, . . . , iV — 1}. We consider 
two such measures. The first measure pi is uniformly distributed on the set B^v, i.e., 



Pl (f) = 2- N \/feM N . 

The second measure P2 is uniformly distributed on the set of results, i.e., 

P2(/) = ajWTT) if a ' = T 

For the average-probabilistic error we want to estimate 



e* vg - pro (M,p) = V pJf) min max \a f - a f (j)\ for i = 1, 2. 

^— ' A: n(A,f)>p j£A 



For the measures p^, the mean of the random variable a/ is clearly |. However, their first 
(central) moments are very different. As we shall see, the moment for the measure pi is small 
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since it is of order iV -1 / 2 whereas the moment for measure P2 is roughly |. Since the first 
moments are the same as the error of the constant algorithm cif(j) = |, we can achieve small 
error of order iV _1//2 for the measure pi without any quantum queries, while this property is 
not true for the measure P2. 

We now consider the measure pi. It is interesting to ask if the QS algorithm has the 
same property as the constant algorithm. We shall prove that this is indeed the case iff M 
is divisible by 4. 

We compute the first moment or the error of the constant algorithm, which is 

TV 

k=0 

We do it only for odd N since the case of even N is analogous. We have 

N 





1 


k 




2 ~ 


N 



E 2 ~ 

k=0 



N 



k 

N 



[N/2} 

2 E 

k=0 



2 n-i_ 2 _ 



Thus 



N 



D avg-pro ( ' 
'Pi 



'(0,1) = ^2 

fc=0 

By Stirling's formula 

k\ = \Phxk 



-N 



1 

2 ~ 

)JV-1 



k 
N 



k 
N 



L/V/2J 

= E 

k=0 

N - 1 
(iV-l)/2 



LJV/2J-1 

2 E 

k=0 

N — 1 
(N - l)/2 



iV- 1 



is even. 



9 k /12k 



we estimate the both binomial quantities in 

1 1 



2-%/-^) if N is odd, 
2-(^)(^ 2 ) ifiVie 



for certain 9k € [0, 1] 

by 

3 l/(12(jV-l)) 



(22) 



proving that 



D avg-pro 
'Pi 



(o,i) = E^C) 

fc=0 v 7 



iV 



27r ViV 



1 , ,„ 1 
l + o(l))< 



1 



3 l/(12(jV-l)) 



2tt v 7 ^" 



(23) 



We are ready to analyze the average-probabilistic error of the QS algorithm. 

Theorem 4. Assume that M is divisible by 4 and let p 6 (|jnr]- 27ien £ 
probabilistic error of the QS algorithms with respect to the measure pi satisfies 



^avg-pro,' 
'Pi 



(M,p)<mmt~,J£-\ll + 




7T- 



,l/(12(jV-l)) 



4 M V 2i 



< 7 7r(l + o(l))miiJ — 



4M 2 v/jV^I 
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Proof. The estimate e^ g - pro (M,p) < e^ g - pro (M, < e™ r " pro (M, ^) is obvious and, ap- 
plying Corollary we get 

3 7T 



3 avg-pro 
'Pi 



(M, A) < 



4 M 



As before denote a a = — arcsin y/a. Let a — \ + x. We are interested in the behavior of 
ax for |x| < 2- Clearly ax = \M. Let \x\ < ~, By Taylor's theorem, we have 



x 



M M 



for ^e(i,i + x) 



and 2a/ (1 — £:r)£r > Vl — 4x 2 . Assume additionally that 

M \x\ 



1 

< - 



7T VI - 4x 2 ~ 4 



which is equivalent to assuming that 

\x\ < 



7T 



(16M 2 + 4tt 2 ) 1 /2- 



Since M is divisible by 4 then 
yields 



(71 



r+2 



~M for x > 0, and 



2" 



mm 



471 j. 

2+2 



°1, I 17 !, _ 
2 2 





1 M 




K 







7T VI - 4X 2 



±M for x < 0. This 



(24) 



Observe that 



(71 



^ G[0 ,J]U[|,1]. 



Indeed, for x < we have 



(Tl 



7+2 



M Ixl 



(Tl 



and for x > we have 



(71 

2" 



_ M |x| 



e[i-ii] = [|,i], 



as claimed. 

Let a = i + x. By the proof of Corollary El the error of the QS algorithm satisfies 



7T 

a — a < — min 

1 1 ~ M 



(71 



7+2 



(7i , a i 



(71 



7+2 



and by (J2U), we have 



\a — a\ < 



x 



VI - 4x 2 
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We split the sum that defines e^ g pro (M, Jj) into two sums. The first sum is for / G B^v 
for which a = aj = \ + x with \x\ < 7r/(16M 2 + 47r 2 ) 1//2 and the second sum is for / for which 
a = df = | + x with |x| > 7r/(16M 2 + 47T 2 ) 1 / 2 . For the first sum we estimate the error of the 
QS algorithm by |rr|/\/l — Ax 2 and for the second sum, by the worst-case error 37r/(4M). 



Hence we have 



3 avg— pro 
'PI 



(M, £) < 



E 



pi(/) 



i°/-§i 



/: |a / -i|<7r/(16M 2 +47r 2 ) 1 /2 



l-4(a / -|) 2 



E Pi(/)t£- 



/:|a / -||>7r/(16A/ 2 +47r 2 ) 1 /2 

Since a/ = fc/A^ for some integer e [0, A], 



4M 



O avg-pro 
Pi 



(Af, 4,) < 



E 



-A' 



+ 



k: |fc/Ar-i|<7r/(16A/ 2 +4 7 r 2 ) 1 / 2 
3 7T 



4M 



E 



-TV 



1 - A{k/N - |) 2 

iV 



fc: |fc/Af-i|> 7 r/(16Af 2 +47r 2 ) 1 / 2 

Since 1 — A(k/N — \) 2 > 1 — 7r 2 / (4M 2 ) > |, the first sum can be estimated as 



E 



fc: |fc/Af-i|<7r/(16A// 2 +47r 2 ) 1 / 2 



N 



\k/N-\ 



l-A{k/N-\) 2 
2 



< 



E 



-N 



k: |fe/AT-i|<7r/(16M 2 +47r 2 ) 1 / 2 



N 



1 k 
2~N 



The second sum can be estimated by 



3 7T 

AM 



E 

fc: |fc/AT~||>7r/(16Af 2 +47r 2 ) 1 / 2 



-N 



N\ (16M 2 + 4vr 2 ) 1 / 2 
fey 7r 



k 1 

iV ~ 2 



< 3 A 1 + 



AM 2 



E 



-N , 



k: |fc/Ar-||>7r/(16Af 2 +4vr 2 ) 1 / 2 

Adding the estimates of these two sums we obtain 





k 


1 





N ~ 
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3 avg-pro f 
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/ 2~~ ^ 

< M '^ 3 V 1 + iM5E 2 



fc=0 
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2 ~ 


N 



The last sum is given by (J22|) and estimated by (|2l?jl . Hence 

3~ 



e avg-pro (M) ^ < 

which completes the proof. 




7T 



1 



2ir V 1 + 4M 2 VN- 1 



,1/(12(JV-1)) 
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In the next theorem we consider the case when M is not divisible by 4. 



Theorem 5. Assume that M > 4 is not divisible by 4, and let p G (|, Then the average- 
probabilistic error of the QS algorithm with respect to the measure pi satisfies 



3 avg— pro 
Pi 



(M,p) > 



71 



4MV 1 "^~^ ),1 ~ 2eXP 



Ntt 2 



ip My 



V/3 > 1. 



Proof. Let M = AM' + r for r 6 {1, 2, 3}. Let, as before, a a = (M/ir) arcsin ^/a. As in the 
proof of Theorem for \x\ < \ we have 

MM x /ii \ 

en = — H . for ^ G U, 5 + x) 

2+* 4 vr 2^(1 



and 2-^/ (1 — £ x )£ x > Vl — 4a; 2 . Assume additionally that 



M \x\ 1 



vr VI - 4x 2 ~ 4/?' 



which is equivalent to assuming that 

\x\ < 



7T 



(16M 2 /3 2 + 4tt 2 ) 1 /2 
Thus for x satisfying (|23|) we have 

., r M x6(x) 



(25) 



and 




= M' and 










2 +x 



(71 = M' + - + t . 

2 +x 4 7T Vl - 4X 2 



M' + l. 



with 0(x) G [0,1], 



From the proof of Corollary [T] we have \o~ af ~\ , |_cr a/ J }, /) > Jy. Since f) > p> 
then either | c a .] G A or [cr a J G A. We then estimate 



e av g - pro(M) ^) > ^ pi(/) min { | tt/ _ g/ ( ) | ; | a/ _ a f ( \a af ] ) |}. (26) 



We now estimate the error of the QS algorithm for / G Bjv such that a/ = | + a; for x 
satisfying (|23j) and the outcome j = \cr af \ = M' or j = [c af \ = M' + l. Denote the outcome 
by M' + k for k G {0, 1}. By Taylor's theorem we have 



sin 



M' + k \ . 2 7T / i < 

— — — TV ) — sin — + — (k — tt 
Ml V4 M v 4 ' 



~ + sin(2^ T )^(K - ±r) for G [|tt 3 Jtt + (tt/M)(k - |r)]. 
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Since sin(t) >2t/n for t 6 [0, 7r/2], we have | sin(2£ KjT )| > 1 — \4k — t\/M. Consider the error 
for the outcome M' + k and x satisfying ()25|). Then |x| < 7r/(4/3M) and the error can be 
estimated by 



+ x — sin^ 



M' + K 

M 



■ 7T 





X 
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M 




7T 


> 






4 



7T 



4k - t 






4k - r 




M 





7T 



4/3M' 



Clearly, |4/c-r| G {1,2,3} and \Ak-t\/M e [1/M,3/M]. Then \4k-t\(1-\4k-t\/M)/M > 
(1 - 1/M)/M. Therefore 



— h x — sm 
2 



'M' + K 

M 



■ 7T 



M 14 



M / 4/5 



7T 



4M 



M /3 



(27) 



Hence, for / such that af = | + x with x satisfying (|25|) we have 

min{|a / -a / (L<T 0/ J)|,|a / -a / ([a- 0/ l)|} > ^ - - -J . 
We are now ready to estimate e^ i e ~ pro (M,p). First, by (|26j). we have 
e av g - pro(M)p) > J2 p 1 (/)min{|a / - 5/ ( [a., J )|, \a f - a f ( \a af ])\}. 

f: |a / -i|<(16Af 2 /3 2 +4 7 r 2 ) 1 /2 

This, (EZ|) and the Bernstein inequality, J2 k . lk/N h >E 2 ~ N ffi ^ 2 e - 7Ve2 / 4 , yields 



O avg-pro 
"PI 



> J] pi(/) 

/: |a / -i|<7r(16A/ 2 /3 2 +47r 2 ) 1 /2 



7T 



4M 



1 1 

M ~ P 



TV 



AM 



M (3 



i-JV 



fc: |fc/V-i|<7r(16M 2 /3 2 +4 7 r 2 ) 1 /5 



> 



> 



7T 

4M 

7T 

4M 



1 - — - 4 J ( 1 - 2 exp 



Af /? 
1 1 

M ~ P 



1 . - — - - ] ( 1 - 2 exp 



4(16M 2 /3 2 - 4) 
iVvr 2 



which completes the proof. 



□ 



Obviously, in the average-probabilistic setting, we should use the QS algorithm with M 
divisible by 4. Then Theorem states that the error is of order min {M~ l ,N- 1 / 2 }. Recently, 
Papageorgiou [TU] proved that for any quantum algorithm that uses M quantum queries the 
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error is bounded from below by of cmin{M _1 , N~ 1 / 2 } with probability p G (|, Here, c is 
a positive number independent of M and N. Hence, the QS algorithm enjoys an optimality 
property also in the average-probabilistic setting for the measure pi as long as we use it with 
M divisible by 4. 

We now turn to the measure P2. Clearly, the average-probabilistic error of the QS 
algorithm is bounded by its worst-probabilistic error, which is of order M~ l with probability 
p £ (f> It turns out, again due to a recent result of Papageorgiou [E] that this bound is 
the best possible, since any quantum algorithm that uses M quantum queries must have an 
error proportional at least to M _1 . Hence, the factor iV~ 1//2 that is present for the measure 
Pi does not appear for the measure P2, and the behavior of the QS algorithm is roughly the 
same in the worst- and average-probabilistic settings. 
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